Detecting Cosmic Gravitational-wave Background from Super-heavy Cosmic Strings 

with LISA 
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Although cosmic string scenario for galaxy formation is disfavored by CMB data, it is of great 
interest in the generation of cosmic gravitational-wave background. This research aims to develop 
an algorithm to extract cosmic gravitational-wave background produced by cosmic strings from the 
LISA data stream, and apply the algorithm to the simulated data stream containing the background 
produced by cosmic strings with various strength to study the detection threshold for this source. 
For 1-yr observation, It is found that the detection threshold of Gfi is 3.12 x 10 -16 in the standard 
scenario. In the case that p and e are adjustable, the detectable region in parameter space is defined 
by (Gm) 2/ VV 1/3 > 4.6 x HT 11 . 



PACS numbers: 04.30. Tv, 98.80.Cq 



I. INTRODUCTION 

The cosmic gravitational-wave background (CGB) is 
the superposition of relic gravitational waves generated in 
the very early universe. Several scenarios might produce 
CGB have been introduced fiL including inflation [hjg , 
vacuum-bubbles colliding [Sj-Qjl], and cosmic strings |15r - 
Inflation scenario, favored by the Cosmic Microwave 
Background Radiation (CMBR) data, introduces pertur- 
bations in all fields, so it is expected the CGB generated 
by inflation would exist. However, the strain amplitude 
of CGB given by this scenario is below the LISA sensi- 
tivity. As for vacuum-bubble colliding and cosmic string 
scenarios, they are the product of the phase transition. 
Though both scenarios could create detectable CGB by 
LISA, in this paper we will focus on the latter case. 

Contrary to localized sources such as binaries, the cos- 
mological sources are randomly distributed across the 
sky and the signals are incoherent, causing an unre- 
solved average of the variation of arm-length and pro- 
ducing a continuum which is entangled with instrumental 
noise. Moreover, the existence of cosmic gravitational- 
wave backgrounds is still unclear. Therefore, the algo- 
rithm used to extract CGB is demanded not only being 
able to separate CGB from instrumental noise, but also 
able to decide whether specific kinds of the backgrounds 
appear in a data set. 

The proof mass will be drifting away from its free- 
falling trajectory due to persisting perturbation by ran- 
dom acceleration. Despite that the drift itself cannot 
be modeled because its behavior is analogous to random 
walk, the trend of the drift can be modeled by a quadratic 
function of time at 2 + bt + c where a, b, c are parame- 
ters need to be fitted. Having removed the trend from 
the time series data and then Fourier transformed the 
residual, one can obtain the noise spectrum in the fre- 
quency domain. The low frequency part of the spectrum 
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is regarded as inherent acceleration noise. In case that 
CGB existing in the data, inaccurate estimate might be 
induced in data analysis if the trend and the background 
are sequentially extracted from the data stream. To avoid 
the bias, the analysis shall be performed entirely either 
in the time domain, or in the frequency domain. The dif- 
ficulty to deal with the both simultaneously was that the 
stochastic nature of CGB obstructs the extraction from 
the time domain, but the function of the power spectrum 
of the trend in the frequency domain was unknown. It is 
natural to extract the CGB in the frequency domain be- 
cause the power spectrum with respect to a specific kind 
of background has its own feature, making itself able to 
be distinguished from others. To successfully separate 
the background and the trend in a data set, we have de- 
rived the function of the power spectrum of the trend 
induced from the random acceleration in the frequency 
domain [l8| so that the trend can be estimated simulta- 
neously with the background. 

Data analysis is a time-consuming task. In general, the 
time required for completing estimation with brute force 
is roughly exponential increasing with the number of pa- 
rameter. Suppose one has a model with 10 parameters. 
Even if one just has ten trial values for each parameter, 
10 billion parameter sets are produced in order to find 
the best estimate. To solve this, the Metropolis algo- 
rithm fl9| has been adapted in our algorithm 18], which 
is a kind of Markov Chain Monte Carlo (MCMC) method 
to estimate parameters. The time used for parameter es- 
timation with a MCMC is roughly proportional to the 
number of parameters. Furthermore, we have applied 
the simulated annealing method to speed up the search- 
ing of the Markov Chain [Io|-|22|. The Gelman & Rubin 
method was chosen as a diagnostic for the convergence 
of the chains to confirm the robustness of the estima- 
tion [23|. In our previous work |l8j, the algorithm has 
been constructed based on Bayesian statistics with the 
three techniques to perform parameter estimation. Here, 
we will extend the capability of the algorithm to per- 
form model selection task, and use it to analyze data 
containing CGB from cosmic strings to study the detec- 
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tion threshold of cosmic string scenario for LISA. 

The organization of the paper is as follows: in section 
ITT1 the background of cosmic string scenario is quickly re- 
viewed, including the gravitational radiation from loops 
and stochastic gravitational-wave background from cos- 
mic strings. In section IIIII we introduce the algorithm 
for model selection. In section HVl we utilize the model- 
selection algorithm to study the detection threshold of 
cosmic string scenario for LISA, and the analysis results 
will be presented. 



II. SUPER-HEAVY COSMIC STRINGS 



the reconnection probability of intersecting superstrings 
may be much smaller than 1, while the ordinary strings 
always reconnect (p=l). The two concerns lead to the re- 
examination of g ravitational radiation from these cosmic 
superstrings jl7j |. 



B. Gravitational Radiation from Loops 

Suppose the trajectory of the string is x(ct, t), a is a 
parameter along the string. The equation of motion for 
a string [l5[ is 



In this section we will quickly summarize previous 
works, mainly by T. Vachaspati and A. Vilenkin, on 
cosmic string [ill, [TtJ ■ The expression for stochastic 
gravitational-wave background from cosmic strings used 
in this paper will be given. 



= 
= 
1. 



(1) 



Dots and primes stand for derivatives with respect to t 
and a. A simple solution with only two frequencies is 



A. Background Review 

Cosmic strings are topological defects which produced 
during a symmetry breaking phase transition in the very 
early universe. Although cosmic string scenario has been 
disfavored by CMBR data, it could be a significant source 
of gravitational-wave background. All strings are in the 
form of loops. In a Hubble volume, there are a number 
of small loops, while there are few loops larger than a 
Hubble volume so that only a section of loops contained 
in the volume, which is called infinite strings. An ex- 
ceptional feather of cosmic strings is that the number 
of infinite strings roughly remains constant during the 
evolution of the universe. The string networks evolve 
as the universe expands These loops are characterized 
by mass density fi. In particular, the tension of loops 
are equal to /i, which is much larger than the tension of 
daily-use strings. Under such strong tension, the loops 
oscillate in a relativistic speed, which leads violent varia- 
tions of mass distributions with respect to time, produc- 
ing gravitational waves with typical frequency / ~ L —1 
where L denotes as their lengths. The cosmic string net- 
works develop with the expansion of the universe. As 
the consequence, the number density of infinite strings 
roughly keep the same in a Hubble volume. Small loops 
would disappear due to releasing energy by gravitational 
waves. However, more loops would be produced from re- 
connection mechanism to replace the absent loops. T. 
Vachaspati and A. Vilenkin have made a careful calcu- 
lation of gravitational radiation from oscillating closed- 
strings [15| where the contributions from each loop in a 
wide range of frequencies was considered. 

Subsequently, it has been suggested that in brane infla- 
tion theory the fundamental (F-) and D-string networks 
could be produced during the evolution of universe. The 
first new character of the superstrings is that their typical 
size may different from that of ordinary strings. Secondly, 
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where plus and minus represent the left and right-moving 
modes of oscillations on loops. This gives a fair represen- 
tation of closed strings produced during the early uni- 
verse. As high-frequency modes are decayed during the 
expansion of universe, the cutting off higher-order terms 
in Fourier series is a fair approximation. The distribution 
of gravitational radiation generated by a closed string is 

^ = ^ (V;>„,fc)T^>„,fc) - i|T>„,fc)l 2 

(3) 

where u> n — Aim/L is the angular frequency, and 
T)n/(wn,fc) is the Fourier transform of the stress-energy 
tensor of a string 

V(x, t) = M J da (x^x y - x V) (x - x(a, t)), 

(4) 

\x is mass per unit length. The total power of the gravi- 
tational radiation produced by the string is 



where 



(5) 



(6) 



It can also be expressed by 7G/Z 2 , where 7 ~ 50 is derived 
from numerical computation. 



C. Stochastic Gravitational- wave Background from 
Cosmic Strings 

The stochastic gravitational background was produced 
by loops with different shapes and sizes around us. The 
formula of total gravitational radiation by a loop has been 
derived. The energy density of the gravitational back- 
ground can be calculated by integrating the product of 
the total power from a loop and the number density over 
entire history of universe, as long as the number density 
of loops is suggested. 

Loops only emit gravitational waves at specific fre- 
quencies f n = 2n/L(t) relating with their length. The 
L(t) = L n — jGfit is their length, where "fG/j, represents 
their contraction caused by the gravitational radiation. 
Their formation rate dn/dt at time t and their initial 
length L n are proportional to j3t~ A and at, respectively, 
where a and (3 are numerical coefficients. By the rela- 
tion mentioned above, it is shown that the number den- 
sity of strings corresponding to the gravitational waves' 
frequency / observed presently is 
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dn n (t) = a 3 /3 



a(L n /a) 
a(t) 



dL n , . 

Li ' {<) 



where the scale factor a(t) represents the expansion of 
universe. The present energy density of the radiation 
with frequency / is 



dE 



to 



a(t ) 



dt dn r 



(8) 



L„/a 



The lower bound of the integration is the time t 
when the loop was produced. 

Wyman, Pogosian, and Wasserman found the con- 
straints from WMAP and SPSS data on the fraction of 
cosmological fluctuations by cosmic strings [2J]. Their 
results suggested that cosmic strings can account for up 
to 7 (14%) of the total power of the microwave anisotropy 
at 68 (95%) confidence level. The corresponding bound 
on G\i is 3.4 (5) x 1(T 7 at 68 (95%) confidence level. M. 
V. Sazhin, O. S. Khovanskaya, M. Capaccioli, G. Longo, 
J. M. Alcala, R. Silvotti, and M. V. Pavlov found two 
nearly identical galaxies are having angular separation of 
1.9 arc sec, which suggests G/i ~ 4 x ICG 7 if the gravi- 
tational lensing is caused by a cosmic string [25| . Janet 



et al. [26| developed a technique to detect a stochas- 



tic gravitational- wave background by finding correlations 
between pulsar observations. Based on this work, a 
method to place an upper bound on the power of a spe- 
cific stochastic gravitational-wave background (according 
to the corresponding model function) by using observa- 
tions of multiple pulsars was presented [13]. Eight year 
of millisecond pulsar timing observations has provided a 
limit on the critical density of gravitational wave from 
cosmic strings: VLawh 2 < 1.9 x 10~ 8 . 

The stochastic gravitational-wave background from 
cosmic strings is composed of the waves emitted by the 
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FIG. 1: The black (solid) curve shows LISA sensitivity. The 
red (dot) curve represents the confusion noise background. 
The green (dash), blue (dash-dot), yellow (double-dot dash), 
and purple (double-dash dot) lines illustrate the the ampli- 
tude for G/i of 10" 7 , 10~ 9 , 10" 11 , and KT 13 , respectively, 
in super-heavy cosmic strings scenario, e and p are kept as 
1. The figure suggests that LISA may be able to detect the 
gravitational waves in the scenario as Gfi is down to 10 -13 . 



string loops across the sky. In [L7j the stochastic back- 
ground is referred to 'confusion noise background' which 
is defined by 



h 2 f ■ 

confusion 
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dz 



n(f,z)h 2 (f,z)Q(n(f,z)-l). (9) 



n(f, z), the number of the string loops in the universe at 
frequency / and redshift z, is given by 

n(f, z) = 100-(/t )- 5/3 a- 8/ Vn(^)G(z), (10) 
P 

and the dimensionless amplitude of gravitational wave 
produced by a single loop is 

h(f, z) = Gva 2 / 3 (ft )- 1/3 MzW - 9 m [a, f, z\) (11) 

where a parametrises the typical size of a loop, to is 
present cosmological time 4.42 x 10 17 second, and 
z eq = 10 3 is the redshift of equal matter and radia- 
tion densities. Here 

<p n (z) = z 3 (l + z)-V°(l + — )"/ 6 , (12) 



(Z) = Z- 1 (1 + Z)- 1 /3 (1 + ^ ) -l/3 ) (13) 



and 



C(z) = 1 + 



9z 



(14) 



is an interpolation between matter-dominant era and 
radiation-dominant era. Bear in mind that the step func- 
tion 0(1 — 9 m [a, /, z]) serves as a cut-off to the Fourier 
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FIG. 2: The black (solid) curve shows LISA sensitivity. The 
red (dot) curve represents the confusion noise background. 
The green (dash), blue (dash-dot), yellow (double-dot dash), 
and purple (double-dash dot) lines illustrate the the ampli- 
tude for e of 1, 0.1, 0.01, and 0.001, respectively, in super- 
heavy cosmic strings scenario. G/i is kept as 10~ 7 , and p 
is fixed as 1. The figure displays that e does not affect on 
gravitational-wave amplitude very much. 



FIG. 3: The black (solid) curve shows LISA sensitivity. The 
red (dot) curve represents the confusion noise background. 
The green (dash), blue (dash-dot), yellow (double-dot dash), 
and purple (double-dash dot) lines illustrate the the ampli- 
tude for reconnection probability p of 1, 0.1, 0.01, and 0.001, 
respectively, in super-heavy cosmic strings scenario. G/i is 
kept as 10~ 7 , and e is fixed as 1. 



component whose mode number is less than one in which 
6 m [a,f,z] is given by 



t (a, /, z) = (aA,r 1/3 (l + ^) 1/6 (1 + — ) 1/6 - (15) 

Zf=n 



The step function Q(n(f,z) — 1) constrains the inte- 
gration over z in the case of multiple sources, which are 
regarded as the origin of confusion. However, this con- 
sideration is not appropriate for LISA. This is because 
LISA is designed to detect continuous sources. In order 
to reconstruct short-duration events such as bursts by 
triangulation, the minimum number of independent de- 
tectors is 3. For the reason, we adopt root-mean-square 
definition 

r a z 

Pstringtf) = h 2 rms (f) = J — n(/, z)h 2 (f, z), (16) 

averaged over all events, to calculate the power spectral 
density of the stochastic background. Substituting Eq. 
([10]) and CU} into Eq. ([16]). we have 
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where to = 4.4 x 10 17 and c = 1. The integration can be 
worked out numerically, and its value is 0.21. Then the 
Eq. (fl~7)) can be reduced to 



PstringU) = 3.65 X 10 



-12 



(Gm) 2/3 

pcV3 



/ 



-7/3 



(18) 



Figure [TJ [5J and [3] illustrate the influences of the three 
parameters, Gfj,, a, and p, on the strain amplitude h. 



III. ALGORITHM FOR MODEL SELECTION 

Suppose that there are two models Mi and M 2 which 
are the possible explanations for a data set. The ratio 
P(M 1 \D,I)/P(M 2 \D,I), namely the Bayes factor B 12 , 
is regarded as the relative preference of Mi over M 2 . 
Determining the better explanation is identical to find 
the Bayes factor. 

Suppose that the two models M\ and M 2 involve r pa- 
rameters pi, ■ ■ ■ ,p r and s parameters qi, ■ ■ ■ ,q s , respec- 
tively where s > r. With the help of Baye's theorem, 
the relative probability P(Mi\D, I)/P(M 2 \D, I) can be 
decomposed to the product of relative likelihood function 
and prior 

_ P(Mi\ D,I) _ P(D\ Mi, I) P(M 2 \ I) 
12 P(M 2 \D,I) P(D\Mi,I) P(M 2 \IY { ] 

The P(M\\I) is the prior of model Mi being true for 
given background information. Here the evidence is can- 
celled out because the identical data set is used. Since 
we do not prefer any specific model before we analyse 
data, the prior ratio P(M 1 \I)/P(M 2 \I) can be set to 1. 
P(D\ Mi, I) is marginal likelihood function 

P(D\ Mi, J) = J P{D,pi,--- ,p r \Mx,I) d Pl ---d Pr . 

(20) 

With Baye's theorem, P(D,pi, ■ ■ ■ ,p r \M\,I) can be fur- 
ther separated to the product of likelihood function 
P(D\ pi,-- - ,p r ,Mi,I) and prior P(pi,--- ,p r \ Mi, I). 
Hence, Eq. (l20l) can be expressed as 



P(D\Mi,I) = J P(D\ Pl ,--- ,Pr,M l3 I) 

X P{pi,--- , Pr \ Mi, I) dpi- --dp,.. (21) 
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Presumably, any particular set of values for 
(pi,--- ,p r ) is not favoured before analysis. Thus, 
the prior can be assigned by a very wide multiple 
Gaussian distribution 



P(pi,--- ,Pr| M it I) 



1 



1 



exp 
exp 



j__eL\ 

X 2S 2 ) 

{-£}• 



(22) 



The likelihood function P(D\ pi, ■ ■ ■ ,p r , Mi, I) can be 
approximated around its best estimate (pi, ■ ■ ■ ,p r ) 



P{D | px,... , Pr ,M r ,I) ttP(D\p- u 



1 



x exp^--(p-p) ax (p-p) 



,Pr,Mx,I) 

(23) 



where <j\ is the covariant matrix, p = (pi, ■ ■ • ,p r ), and 
p = (pi,- - ■ ,p~r)- Since the uncertainty of the prior is 
very board, the centre of prior can be shifted to the best 
estimate 



1 1 

P(px,--- ,p r | Mx,I) « X ••• X 



2tt(5i 



27T(5 r 



exp{-i(p-pf A-^p-p)}^) 



Here the A -1 is a diagonal matrix where the z-th term 
is Sf. Multiplying the prior and the likelihood function, 
we can get 



P(pi,--- , Pr \ Mx,I)xP(D\ Pl ,--- ,Pr,M X ,I) « P(^|pi,-- - ,^,Mi,I)x 



2tt5i 



27T(5 r 



exp { - i(p - pf K 1 + 2A" 1 )(p - p)}. 



(25) 



A 1 only affects the diagonal term of a , and the i-th since Su ^> <j\ a. Hence, Eq. (|2"5|) can be approximated 
diagonal term of the matrix erf 1 + 2A~* is to 



2 _Sf i + 2o* 



r 2 S 2 a 2 ft 2 

'l it °ii °l 



(26) 



P(pi,--- ,Pr\ Mi, I) x P(D\ pi,-- - , Pr ,Mi,I) a P(D| pi,-- - ,p r ,Mi,I) X -=- 

V 27T0l 

x exp { ^ p) Tfj r 1 ( p ~ p)}- 

I 



27T(5 r 



(27) 



To evaluate the integral in Eq. (f2"Tj) we need to use the 
integral of r-dimensional multivariate Gaussian: 

J ■■■ J cx p { - \{p - p) T o-r x (p - p)} dpi ■ • ■ dPr 

= (47r) r /V^ei(<7i). (28) 



With Eq. {28]), substituting Eq. (J27J) into Eq. (J2TJ) we 
can obtain 



P(fl|Jfl ,, )s ^ f(fl | fl , 



,Pr,Mi,I). 



Sx • • • S r 

(29) 

In principle, the values of Sx,-- - ,S r are considered as 
infinite to satisfy the condition of the prior being very 
board. In practice, there is a way to assign fairly large 
values for the uncertainties of prior. For instance, the 



uncertainties can be set by the range of the parameters, 
which can be known either from the experiment design 
or the results of previous experiment, or guessed from 
theoretical calculation. 

By repeating the same steps, we can have the expres- 
sion for P(D\ M 2 , I), and then we can acquire the relative 
probability of M% and Mi 



P{D\ M 2 ,I) 
P(D\ Mi, I) 



\ / 2 s - r 



Si--- 8 r \/Det{a 2 ) 



S'i 



P(D\pi 



5' s ^Det{ai) 
,Ps',M 2 ,I) 



P(D\p U --- ,p r ,Mi,I) 



(30) 



where the symbol prime refers to the objects related to 
M 2 model, and Ci and a 2 are the covariant matrixes given 
by fitting with Mi and M 2 model respectively. Since the 
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TABLE I: Bayes factor confidence level cited from [2S 



as 



B12 


Evidence for model Mi 


< 1 


Negative 


1 ~ 3 


Not worth more than a bare mention 


3 ~ 12 


Positive 


12 ~ 150 


Strong 


> 150 


Very Strong 



noise is independent and additive, the likelihood func- 
tion P{D\px, ■■■ ,p s ', M 2 , 1) and P{D\pi , • • • ,p r , M u I) 
can be decomposed into a series product of Gaussian dis- 
tributions. Table IIIII lists the evidence for model Mi 
against the Bayes factor. 

If the noise a is unknown, we can integrate the series 
product of Gaussian distributions over a to obtain the 
marginalised likelihood function, which is 



N 



P{D\p u --- ,p r , Mi, I) 

(31) 

where F(pi,---,p r ) is defined as J2i=i(Pi ~~ 
p{ {pi, ■ ■ • ,p r ))- The other marginalised likelihood 
function P{D\pi , ••• ,p s , M2,I) is given by the same 
step. With the two marginalised likelihood functions, 
we have the likelihood ratio 



P(D\pi 



■P. 



;,M 2 ,I) ( F(pi,--- ,p r ) \ N ~ 1/2 



P(D\p u --- ,f r ,Mi,I) \F(fi',--- ,p s ')_ 

(32) 

in the case that the noise is unknown. If the correla- 
tions among parameters are small, U Det(p\) can be fur- 
ther approximated to successive product of uncertainties 
0pi • • • o~ Pr . We define the small as the magnitude of cor- 
relation being smaller than 0.2, in which case the dis- 
crepancy between the approximation and the true value 
of determinant is about 5%. Expanding \J Det(a 2 ) as 
cr x • • • x a ps , the first term on the right hand side of 
Eq. (|30|) can be simplified as 



P(D\ M 2 ,J) _ r- v Pr+1 
P(D\M U I)~ V - 



■V+l 



P(D\Pi,--- ,p' r+ i,M 2l I) 
P(D\pi,--- ,p- r ,Mi,I) ■ 



(N-l) \fF(pi,--- , p - r )^-(N-i/2) 
N N 



. y /^~T 5i---5r\/Det{o 2 ) = jh Sr_ 

5'i---S' s ^Det{ai) <J P1 ' " 



(33) 

This term is the so-called Ockham factor. To understand 
this factor, we suppose that the model M 2 is built upon 
Mi with an extra parameter p r +i- The uncertainties of 
parameters estimated by different models are roughly the 
same, so the right hand side of Eq. (f3"3"|) can be approxi- 
mated by 



— 5i---5 r ^Det{o 2 ) 



(34) 



S' 1 ---5' sy /Det(a 1 ) 
The relative probability of M 2 and Mi can be simplified 



(35) 

The second term in Eq. (|35p is the likelihood ratio of two 
models. Since a more complex model will yield a better 
fit than a simpler model, the likelihood ratio in Eq. ([33)1 
is always larger than 1. On the other hand, it is found 
that the Ockham factor is much smaller than 1 as long 
as the quality of data is not poor, (poor-data limit can 
be expressed as u r+1 3> <5 r +i-) It plays a role of 'penalty' 
for using an extra parameter to fit data. As we can see, 
if the more complex model cannot produce a much bet- 
ter fitting, giving a high likelihood ratio to overcome the 
penalty, the simpler model will be suggested. Although 
a more complex model always fits data better, a simpler 
one might be preferred if the discrepancy of its descrip- 
tions to data is not much worse than the more complex 
ones. 

To make a judgement on the number of signals in data, 
we can use the simplest model Mi which assumes that 
there is no signal contained in data, and the model M 2 
which assumes that there is one signal in the data to 
analyse the data, respectively. Then we compare their 
relative probability P(M 2 \D,I)/P(M 1 \D,I). If the ra- 
tio is smaller than 1, it indicates that the data are just 
purely noisy; otherwise, it is suggested that there may 
be one or more signals out there. We further move 
to the more complex model M3 and compute the ratio 
P(M 3 \D,I)/P(M 2 \D 7 I). If the ratio is smaller than 1, 
it is suggested that there is one signal in the data; oth- 
erwise, we move forward to the more complex model M4 
and follow the same steps to compute the ratio and so 
on until it is smaller than 1. The likeliest model would 
advise us the best description of the data. 

In the end, we would like to discuss the uncertainty 
of prior 5 r+1 in Eq. (1351) deeper. A concern is that the 
uncertainty seems to be free to choose. At least, vari- 
ous approaches would propose different values of the un- 
certainty. Therefore, Mi model may be inclined with 
the choices of some values, but M 2 model may be indi- 
cated with the others. In other words, the model selec- 
tion is worried to be subjective. To clear this unease, we 
should keep in mind that choosing a bigger value of the 
uncertainty just represents the ground is more conser- 
vative. The bigger uncertainty might lead to reject the 
signal which is indicated by a model with a smaller un- 
certainty. However, it is just shown that the signal is not 
strong enough to be accepted by the conservative point of 
view. The relative probability P(M 2 \D,I)/P(Mi\D,I) 
just can be considered as the relative degree of belief 
between the two models. After all, we choose to regard 
probability as the degree of belief rather than a long term 
occurrence rate as we use Bayesian statistics in the first 
place. 
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FIG. 4: On the top plot the black solid curve is the simulation 
of random acceleration and shot noise on the 'sky plane'. The 
red dash line is the cosmic gravitatinal-wave background pro- 
duced by cosmic strings. On the bottom plot the black solid 
curve is the full spectrum including both the background, the 
random acceleration, and the instrumental noise plotted on 
the 'sky plane'. 



IV. DATA ANALYSIS 

To apply the model selection algorithm to determine 
the existence of the gravitational-wave background gen- 
erated by the cosmic strings in a spectrum, we need 
two models to describe the spectrum in which the back- 
ground is present or absent, respectively. Models were 
constructed in the hierarchical sense, which means that 
all models are built on a primary one. In the case of the 
background from cosmic strings being absent, the spec- 
trum only contains the power from the drift trend due to 
random acceleration. Because the trend definitely exist, 
it is chosen as our primary model. The model, namely 
Mo, can be expressed by Eq. (|36]l: 



M (/„) = 



4i 2 



1 



N 2 



1 



sin 4 2™ 3 sin 2 ft 



(36) 



where N, the total number of data, is set as 8192 in our 
research, n is the index. A, the sampling rate, is set 
as 1.5 sec, and /„ = n/NA is the frequency. The total 
observation time NA is 12288 sec. The data from this 
setup can cover the LISA frequency band 0.1 Hz to 0.1 
mHz. L is the arm-length of LISA, and a is the averaged 
acceleration need to be estimated. 

The complex model is established by adding a specific 
kind of signal to the simpler one. In the case of the 
background from cosmic strings being present, the spec- 
trum is contributed from the drift trend due to random 
acceleration as well as the background. As for the back- 
ground, the physical parameters p, e, and G/x in the Eq. 
([T^| only change the strength of the spectrum but do 
not influence the spectrum shape. We cannot determine 
these three parameters uniquely through the strength of 
gravitational- wave background by the cosmic strings un- 
less alternative information is acquired. Hence, we define 



a new parameter A in Eq. (|17p to describe their combined 
influences on amplitude as 

A = 3.65 x io-42 (g/ff Hz 7/ 3 _ (37) 
pe 1 ' J 

Then the expression for analysing the spectrum is re- 
duced to 



string 



(f n ) = a f-y\ 



(38) 



The model for the power spectra including both of ran- 
dom acceleration and cosmic string background, Mi(/ n ), 
can be expressed by 



Mi(/„) 



NA 4 a 2 
4L 2 



1 



N 2 



1 



sin 



N 



3 sin 



N 



1 



R(fn 



7/3 



(39) 

The spectrum produced by the cosmic strings was syn- 
thesised with LISA instrumental noise on the 'sky plane' 
to produce full spectra containing both of signal and noise 
as follows: 



Pd(f) = Patringif) + 



PlrJJ) 

R{f) 



(40) 



where Pins(f) is the simulated power spectrum due to 
the LISA shot noise and the random acceleration, and 
R(f) is the response function of LISA to stochastic back- 
grounds. The Pd(f) is the spectra we used to test the 
algorithm. An example of the simulated full spectrum is 
shown in the Figure HI 

The data analysis starts from parameter estimation by 
the primary model Mo. Then the model M\ is utilized to 
repeat the process. Their results are compared through 
eq. ([55]) . If the model Mo is preferred, it is indicated 
that there is no signal in the data set; otherwise, the 
background from cosmic strings may exist. 



A. Upper Bound for the combined amplitude A 

In the model Mo(/), the magnitude of random acceler- 
ation a is the only parameter. As for Mi(/), the model 
includes two parameters, the magnitude of random ac- 
celeration a and the strength of the gravitational-wave 
background by cosmic strings A, To apply Eq. (|3"0")) to 
compare the model Mq(/) and Mi(/), we have to know 
da/da^A — 1/Sa where 8a is regarded as the allowable 
range for the parameter A , Max(A)-Min(A). Since the 
minimal value of amplitude is zero, to determine 5 A, the 
rest thing we must know is the upper bound of A. 

The parameter A is associated with the energy density 
of the gravitational-wave background through Eq. (4.7) 
in [l7|: 



3tt 2 



(fto) 2 h 2 (f) = 2.87 x lO 36 ^/- 1 / 3 . (41) 




FIG. 5: The green (solid) curve is combined spectrum includ- 
ing both cosmic string signals and instrumental noise. The 
red (dash-dot) curve is the residuals after fitting. The blue 
(double-dot dash) line is the background produced by cosmic 
strings in the spectrum. The black (dash) line is the recovery 
of acceleration noise suggested by fitting with the Mo model. 
The observation duration is 12288 seconds. 



Therefore, giving an upper limit U of fl g (f), the upper 
bound of A can be obtained as 



A < 3.49 x lQ- 37 Uf 1/3 Hz 7/3 . 



(42) 



The upper limit Gy, < 6.1 X 1(T 7 given by the CMB data 
0H1 yields the upper bound A < 2.63 x 1(T 46 Hz 7 / 3 
through Eq. (f37|) . The analysis of 8 year of millisec- 
ond pulsar-timing data provided VL g h 2 < 6 x 10~ 8 at 
/ ~ 1/(7 yr) according to the original analysis (29|, or 
Qgh 2 < 9.3 x 10 -8 based on the Bayesian approach (30j 
where h 2 0.65 2 = 0.42. The two constraints lead to 
A < 8.25 x 10~ 47 Hz 7 / 3 for the original analysis, or 
A < 1.28 x 10~ 46 Hz 7 ! 3 for the Bayesian approach. From 
the LIGO S4 science run the limit on the isotropic back- 
ground for the scale-invariant case is fi s = 1.20 x 10~ 4 
at / = 100 Hz, and the limit for the constant strain 
power case is tt g = 5.13 x 10~ 5 at / = 100 Hz [H|. 



From the upper limits we can obtain the upper bounds 
A < 1.94 x 10~ 40 Hz 7 ' 3 and A < 8.31 x 10~ 41 Hz 7 / 3 
for the two cases, respectively. The result of LIGO sci- 
ence run S5 further constrains the upper limit of fre- 
quency independent O s below 6.9 x 10~ 6 at / = 100 Hz 
[321 ] . giving the upper bound 1.12 x 10~ 41 on A. From 
these constraints we know that the lowest upper bound 
is A < 8.25 x 10~ 47 Hz 7 / 3 , but to prevent the spurious 
background being selected, the more conservative upper 
bound A < 1.28 x 10~ 46 Hz 7 / 3 is chosen for the Max(A). 
Therefore the S A is set as 1.28 x 10~ 46 Hz 7 / 3 . 



B. Results 

Then the existence of the background produced by cos- 
mic strings in the data sets is suggested by the model 
selection method employing the Laplace approximation 



FIG. 6: The green (solid) curve is the full spectrum, and 
the blue (double-dot dash) line is the background fed into the 
spectrum, which are the same with those in Fig [5] The black 
(dash) curve shows the recovered power caused by the random 
acceleration. The yellow (double-dash dot) line represents the 
recovered background produced by cosmic strings, which is 
slightly under the real one. The red (dash-dot) curve is the 
noise part of the data, representing the difference between the 
data and the Mi model. 



as described in Sec. IIIII From model selection mothed, 
the probability ratio of Mi to Mq is given by 



P(Mi\I) = 1 y/Detfa) L(x! \M U I) 
P(M \I) SA y/Detto) X L(xa \M ,I)' 



(43) 



where L is likelihood function, I is background informa- 
tion. Xi and Oi are the best estimate and the covariant 
matrix with respect to model M.%. SA is the bounded 
range of parameter A. Mo and Mi, described by Eq. 
([36]) and Eq. (|39[) respectively, are the model for the cos- 
mic background being absent or present in the spectrum. 
The likelihood ratio L(x\ \Mi, I)/ L(xq \Mq,I) is com- 
puted through Eq. (|3"2"|) . The two models are applied to 
analyse same data set subsequently. The best estimates 
and the covariant matrix with respect to each model will 
be given. Next, the best estimates for the two models 
are used to compute their likelihood ratio. Finally, com- 
bining the likelihood ratio with the determinants of the 
covariant matrices, Eq. (|4*5)l provides the Bayes factor 
of two models. If the ratio is higher than 1 and reach 
to certain level, it indicates that the instrumental noise 
with cosmic string model is the better description of data, 
which is the evidence of cosmic strings existence; other- 
wise, pure instrumental noise will be preferred and the 
existence of cosmic strings will be declined. The ratio of 
20, which means the existence of cosmic strings is 95% 
likelihood, is our definition of LISA detection threshold. 

To pin down the LISA detectability of cosmic string 
spectrum strength, firstly we apply the model selection 
algorithm in the analysis of a full spectrum containing 
both instrumental noise and a background generated by 
cosmic strings corresponding to a specific strength. If 
the Mi model is selected, it means that the strength of 
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TABLE II: Analysis Result of Gravitational-wave Background by Cosmic Strings. The first column A is the real 
value of the background. The second column a is the magnitude of random acceleration fitted by the Mo model. The third 
column a is the magnitude of random acceleration fitted by the Mi model. A in the forth column is the estimation of the 
background strength by the Mi model. The p a A is the correlation between a and A in the Mi model. R and P are the likelihood 
and probability ratio of Mi to Mo , respectively. 





Mo Model 




Mi Model 








A (xl(T 49 Hz 7 ' 3 ) 


a (xlO -15 m/s 2 ) 


a (xlO -15 m/s 2 ) 


A (xlO" 49 Hz 7 ' 3 ) 


PaA 


R 


P 


4.2 


2.213 ±0.006 


2.200 ± 0.007 


2.73 ± 0.56 


-0.48 


4.81 x 10 4 


20.60 


4.5 


2.214 ±0.006 


2.200 ± 0.007 


3.02 ± 0.56 


-0.485 


3.28 x 10 5 


143.26 


4.8 


2.215 ±0.006 


2.200 ± 0.007 


3.35 ±0.55 


-0.48 


2.53 x 10 6 


1.06 x 10 3 


5.1 


2.216 ±0.006 


2.200 ± 0.007 


3.66 ± 0.55 


-0.48 


2.21 x 10 7 


9.28 x 10 3 


5.4 


2.217 ±0.007 


2.200 ± 0.007 


3.98 ±0.55 


-0.48 


2.30 x 10 8 


9.60 x 10 4 



le-19 




Frequency (Hz) 



FIG. 7: The green (solid) curve is combined spectrum includ- 
ing both cosmic string signals and instrumental noise. The 
red (dash-dot) one is the residuals after fitting. The blue 
(double-dot dash) line is the threshold for cosmic string de- 
tection. The black (dash) and yellow (double-dash dot) lines 
are the recovery of acceleration noise and string spectrum 
suggested by fitting with the Mi model, respectively. The 
observation duration is 12288 seconds. 

the background is stronger than the detection threshold. 
In this case we decrease the strength of the background 
to generate a new full spectrum, and then re-apply the 
model selection algorithm in the analysis until the Mq 
model is selected. On the other hand, if Mq model is 
selected, it means that the strength of the background is 
weaker than the detection threshold. In the case we en- 
hance the strength of the background to generate a new 
full spectrum, and then re-analyse the data set until the 
Mi model is selected. The strength of the background 
giving the odds ratio P{Mi)/P(M ) of 20 is the detec- 
tion threshold for the gravitational- wave background pro- 
duced by cosmic strings. 

To begin the study of the performance of the algo- 
rithm in the analysis of the gravitational-wave back- 
ground produced by cosmic strings we apply the algo- 
rithm to the full spectrum with the background strength 
A = 7.7 x 10~ 49 . At first the Mq model is used to anal- 



yse the data set. Since it just contains one parameter a, 
we generate 10 chains to test the convergence. From the 
samples we obtain a = 2.22 x 10~ 15 ±6.79x 10~ 18 m/s 2 . 
The analysis result is shown in Fig [5l From the fig- 
ure we can notice that the noise character is different 
from the LISA sensitivity curve below 10 mHz. Next, 
the Mi model is used to analyse the same data set. 
It includes two parameters, a and A, so we generate 
20 chains to test the convergence. From the samples 
we have a = 2.20 x 10~ 15 ± 7.27 x 10~ 18 m/s 2 and 
A = (6.11 ± 0.57) x 10~ 49 Hz 7 ' 3 . The estimate strength 
is deviated from the true value over 2-cr but within 3-er. 
The correlation between the parameters is -0.48. This 
may be due to the similarity between the formulae for 
the acceleration and the background. Because the power 
contained in a frequency bin is fixed, the power accounted 
for the background will be less if that accounted for the 
acceleration is more, and vice versa. The Fig [JJ] illustrates 
the recovered background estimated by the Mi model. 
The likelihood ratio of Mi to Mq is 3.99 x 10 15 . The 
probability ratio is 1.67 x 10 12 , suggesting that the ex- 
istence of the background produced by cosmic strings is 
preferred. 

In this case the gravitational-wave background is suc- 
cessfully recovered. To look for the detection threshold 
for the background we decrease the strength of the back- 
ground in the full spectra, and apply the model selection 
algorithm to compute the probability ratio of Mi to Mq 
model. 

As shown in the Table HH the cosmic string spectrum 
in the standard scenario (p = e = 1) reaches the de- 
tection threshold if A = 4.2 x 10~ 49 . Fig. [7] exhibits 
the LISA detectability limit for cosmic strings by strain 
amplitude. If the amplitude of cosmic string spectrum 
is higher than the blue line, the existence of the back- 
ground is suggested by the algorithm; otherwise, it is 
disfavored. From Eq. (|3"T|) we obtain that the detection 
threshold of Gfi is 3.90 x 10~ u in the standard scenario, 
and the detectable region in parameter space is defined 
by G/i 2/3 /pe 1/3 > 1.15 x 10~ 7 . 
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